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We outhne a statistical theory of turbulence based on the Lagrangian formulation of fluid motion. 
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The strategy of approaching the phenomenon of fully developed turbulence by considering 
the statistics of Lagrangian fluid particles has a long tradition dating back to the early works 
of Taylor Richardson (for an overview, see 0], [|]). Recently, interest in the 
Lagrangian statistics has been renewed by rigorous results on passive scalar advection in 
the Kraichnan model [Q. Additionally, experimental progress has opened the way to gain 
accurate data on the motion of tracer particles which allows one to directly evaluate the 
statistics of the particle acceleration 0, [|] fl^. These experiments will also shed light 



on the physics of relative dispersion of fluid particles ||TT 



A statistical formulation of the problem of turbulence starting from the Lagrangian point 
of view will lead to considerable progress in modeling turbulent flows by Lagrangian pdf 



method |T2|, ||T5[. This method can successfully deal with passive scalar transport. 



turbulent flows involving chemical reactions or combustion |14|. In this approach, which 
origionally dates back to Oboukhov Fokker-Planck equations are used to model the 
statistical behaviour of fluid particles. Although this is an appealing approach, only few 
efforts have been made to relate the Lagrangian turbulence statistics to diffusion processes 
by a direct investigation of the Navier-Stokes equation. An exception is the work of Heppe 
|2l| , who uses a projector- formalism to obtain a generalized diffusion equation for the joint 
velocity-position probability distribution of one particle. 

In the following we shall present a hierarchy of evolution equations for N-point probability 
distributions describing the behaviour of N particles in a turbulent flow. The hierarchy is 
in close analogy to the one presented by Lundgren and Ulinich and Ljubimov ^4] for 
the probability distributions of the Eulerian velocity field (see |^). Additionally we derive 
a functional equation for a suitably defined probability functional, which is the analog of 
Hopf's functional equation ||22| and, therefore, is a concise formulation of the problem of 
turbulence in the Lagrangian formulation. 

Then we shall adress the question whether a type of diffusion process in the sense of 



the random force method advocated by Novikov can approximate the motion of N 
particles. We shall derive a generalized Fokker-Planck equation involving memory terms, 
which determines the evolution of the N-point probability distribution. However, the drift 
and diffusion terms of the generalized Fokker-Planck equation are expressed in terms of 
conditional probabilities of higher order such that the problem remains unclosed. In a 
following paper we shall adress closure approximations. 
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I. FORMAL LAGRANGIAN DESCRIPTION 



In the present section we shall introduce a formal Lagrangian description of fluid flow, 
which will be suitable for formulating evolution equations for statistical quantities. We 
consider the Navier-Stokes equation for an incompressible Eulerian velocity field u(x, t) : 

d 

— u(x, t) + u(x, • Vu(x, t) = - Vp(x, t) + zyAu(x, t) . (1) 
ot 

In order to obtain a closed evolution equation one has to express the pressure as a functional 
of the velocity field. As is well known the pressure is governed by the Poisson equation 

Ap(x, = - V • [u(x, t) ■ Vu(x, t)] . (2) 

In the case of a finite fiuid volume V appropriate boundary conditions have to be formulated. 
Let us consider the case with a prescribed normal component of the pressure gradient. The 
solution of this von Neumann boundary value problem reads 

p(x,t) = i-|^dxG(|x-x'|)V-[u(x',i)-Vu(x',i)] 

+ ^_^^G(x-x')Vp(x',t)-dA (3) 

Here, G(x — x') denotes the Green's function 

AG'(x-x') = -47r(5(x-x') 
n-Vx'G'(x-x') = x'eSV . (4) 

S is the area of the surface SV enclosing the fluid. The boundary condition is incorporated 
into the definition of the Green's function. 

In the case of an infinitely extended fiuid volume the boundary term vanishes and the 
Green's function is given by 

G(x,xO = ^^ • 
Now we can state the Navier-Stokes equation: 

^u(x,t) + u(x,t)-Vu(x,i) (6) 



= - /" (ix'r(x,x') : u(x',t) : u(x', t) + F(x, t) 
+ V y^dx'L(x,x')u(x',i) 



In order to obtain a convenient formulation for the subsequent statistical treatment we have 
introduced the following notation for the pressure and the viscous terms, respectively: 

1 

[r(x,x') : u(x',t) : u(x',t)]„ = G(x - x')u^(x', t)u^(x', t) 

F(x, t) = f VxG'(x - x')n ■ Vj9(x', t)] ■ dA 

Att Jsv 

L(x,x') = Ax5(x-x') . (7) 

The quantity -L(x, x') is a generalized function and is defined in a formal sense. The term 
F(x,t) is due to boundary conditions, c.f. eq. (|^). Note, that we have assumed that the 
normal pressure gradient is prescribed at the boundary. 

Now we turn to a Lagrangian formulation of the equation of fluid motion. To this end 
we consider a Lagrangian path X(t, y) of a fluid particle, which initially was located at 
^(^O)y) = y- The velocity of the particle is given in terms of the Eulerian velocity field 
u(X(t,y),t), whereas the Navier-Stokes equation takes the form 



|x(t,y) = u(X(t,y),t) (8) 
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^^-u(X(t,y),t) = - J dxT(X(t,y),x') : u(x',t) : u(x',t) 

+ F{X{t,y),t) + uJ rfx'L(X(t,y),x')u(x',t) 
For the evaluation of the integrals we perform a coordinate transformation 

x=X(t,yO . (9) 
Due to incompressibility, the Jacobian equals unity: 

aX^ = 1 . (10) 

Now we define the Lagrangian velocity U(t,y) according to 

U(t,y) = u(X(t,y),t) . (11) 

As a result we end up with the following Lagrangian formulation of the basic fluid dynamics 
equation: 

|x(t,y) = U(t,y) (12) 

|u(t,y)) = - J dy'r[X(t,y),X(t,y')] : U(t,y') : U(t,y') 

+ F(X(t,y),t) + z/| f/y'L[X(t,y),X(t,y')]U(t,y') . 



We have obtained a representation of the acceleration of a Lagrangian particle, although this 
expression is formal due to the appearance of the generalized function L(x, x'). Nevertheless, 
we shall find that one recovers meaningful and well defined expressions when one proceeds 
to a statistical formulation. 

II. EULERIAN AND LAGRANGIAN PROBABILITY DISTRIBUTIONS 

The purpose of the present section is to develop a statistical description of the fluid 
motion. The central quantities will be N-point position-velocity probability densities and, 
as the most general quantity, a velocity-position probability functional. 

We start by defining the N-point Lagrangian velocity-position distribution function 



f^{{uj,Xj,yj};t) = 

<5[xi-X(t,yi)]5[ui-U(t,yi)]....5[x;v-X(t,yAr)]5[u;v-U(t,y;v)] > • (13) 



This distribution function allows one to statistically characterize the behaviour of a swarm 
of N fluid particles. The brackets indicate averaging with respect to a suitably defined 
statistical ensemble. Since we have X(t = 0,y) = y the initial condition at t = reads: 



where (yf^({uj,yj}) is the joint probability distribution for the N velocities of the fluid 
particles at initial time t = 0. In the stationary case, g^{{uj, yj}) is the Eulerian probability 
distribution defined below. 

The N-point probabihty distribution for the particle locations, p^({xj , y^ }; t) is obtained 
by integration with respect to the velocities uj and is the central quantity in the theory 
of dispersion of particles in turbulence. Alternatively, we may integrate eq. (|I^) over the 
spatial variables and obtain the pure velocity probability distributions of the particles. 

It is convenient to additionally define the corresponding Eulerian probability distribution 



These probability distributions are obtained from the corresponding Lagrangian (|T^) ones by 
integration with respect to y^, since the fluid flow is incompressible such that the Jacobian 



r({u,,x,,y,};t = 0) = 

5(xi - yi)....5(x7v - yN)9^{{uj,yj}) 



(14) 



(ui,xi; ..; ujv,xjv;t) =< 6{ui - u(xi, t))..5(ujv - u(xAr,t)) > 



(15) 
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(p^) equals unity. Furthermore, the Eulerian probabihty distributions have to fullfill the 
following consistency conditions 

VxJe (ui,xi; ..;u7v,X7v;t) = (16) 

[-Vx'Vu, ■ y C?u'u7^+^(u',x',Ui,Xi; ..;UAr,XAr;t)]x'=Xi i = l,..,N 

These relationships are a consequence of the fact that the spatial dependency of the proba- 
bility distribution enters via the velocity field u(x, t). 

Finally, we mention that it is also straightforward to define mixed Eulerian and La- 
grangian probability distributions. They can be obtained from the Lagrangian probability 
distribution by integration over the variables y^, for which an Eulerian description is per- 
formed: 

-^li ■■' ^mi ym+li --i X^, y„, 

< 6{ui - u(xi, t))..5(u„ - u(x„, t)) 

5[x„+i - X(t,y„+i)](5[u„+i - u(X(t,y„+i),t)].... 

5[x„, - X(t, y„)]5[u„ - u(X(t, y„), t)] > 

= y"dyi...rfym/^(ui,xi,yi; ..;u„,x„,y„;t) . (17) 

Since the Navier-Stokes equation form a classical field theory a complete description 
of the Lagrangian turbulence statistics is contained in the probability density functional 
P[x(y), u(y), t] defined by an average over functional delta distributions: 

P[x(y), u(y), t] =< P[x(y) - X(t, y)]D[n{y) - U(t, y)] > . (18) 

A probability functional for the field X(t, y') is obtained by functional integration with 
respect to the velocity: 

G[x(y),t] = J Du(y)P[x(y),u(y),t] . (19) 
The functional Fouriertransform of the quantity (|18D is the characteristic functional 

Z[r/(y),k(y),t] =< e^/'^yt''(y)-u(*'y)+'^(y>^(y'*)l > . (20) 
This is the Lagrangian analog of the functional introduced by Hopf . 
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III. HIERARCHY OF EVOLUTION EQUATIONS FOR THE N-POINT PROBA- 
BILITY DISTRIBUTION FUNCTIONS 



In the present section we shall obtain an infinite chain of evolution equations for the 
Lagrangian probability distribution functions /^({uj, x^, yj}; t), = 1,2, .... For the sake 
of simplicity we consider a fluid in an inflnitely extended space. Therefore, we do not have 
to take into account boundary terms. The extention to the case of a fluid in a bounded area 
is straightforward. 

We start by considering the one-point pdf /^(u, x, y;f), for which we try to obtain an 
evolution equation by calculating the temporal derivative: 
d 

— /^(ui,xi,yi;t) = 

^ < (5[xi - X(t, yi)]5[ui - u(X(t, y,),t)] > 

= -Vxi- < X(t,yi)5[xi - X(t,yi)]5[ui - u(X(t, yi), t)] > 

-Vu,- <u(X(t,yi),t)5[xi-X(f,yi)](^[ui-u(X(t,yi),t)] > . (21) 

For the flrst term we can use the relationship 

< X(t,yi)(5[xi - X(t,yi)]5[ui - u(X(t, yi), t)] >= 

< u(X(t,yi),t)5[xi - X(t,yi)](5[ui - u(X(t, yi), t)] >= 

ui < 5[xi -X(t,yi)]5[ui -u(X(t,yi),t)] > . (22) 

In order to evaluate the second term we have to insert the Lagrangian formulation of the 
Navier-Stokes equation (|T2|): 

< u(X(t, yi), t)(5[xi - X(t, yi)]5[ui - u(X(t, yi), t)] >= 
-Jdy'< r[X(t,y),X(t,y')] : u(X(t,y'),t) : u(X(t,y'),t) 
5[xi - X(t,yi)]5[ui - u(X(t,yi),t)] > 
+iy J dy' < L(X(t,y),X(t,y'))u(X(t,y'),t) 

5[xi - X(t, yi)](5[ui - u(X(t, yi), t)] > . (23) 

The aim is to relate this term in some way to a Lagrangian probability distribution. This 
can be achieved by inserting the identity 

1 = Jdu^J rfx2(5[x2 - X{t,y')]S[u2 - u(X(t,yO,t)] . (24) 
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As a result one obtains: 

< u(X(t,yi),t)J[xi - X(t,yi)]5[ui - u{X{t,y,),t)] >= 

- j du2 J (ix2 j (iy2r(xi,X2) : U2 : U2/^(ui, Xi, yi; U2, X2, y2; 

+u j du2 j d-x.2 j (iy2L(xi,X2)u2/^(ui,xi,yi;u2,X2,y2;t) • (25) 

We now combine the above formulas to obtain the evolution equation for the Lagrangian 
probability distribution /^({ui, Xi, yi}; t): 

^ (ui , xi , yi ; + Ul • Vx (ui , xi , yi ; (26) 
= Vui • du2 j dx2 J dy2r(xi,X2) : U2 : U2/^(ui, Xi, yi; U2, X2, y2; t) 
-'^Vui ■ j du2 j dx2 j c?y2L(xi,X2)u2/^(ui,Xi,yi;u2,X2,y2;t) 

Due to the nonlocality of the pressure and the viscous term, which relates the Lagrangian 
path of a particle under consideration to the paths of different particles, the evolution 
equation of the one-particle pdf is linked to the two-particle pdf. 

It is straightforward to prove that the N-point distribution function fuUfills the evolution 
equation: 

|r ({u„ X, , y,}; t) + E u» ■ Vx.r ({u„ x„ y,}; t) (27) 
^J2^^i-J(^^J(^^j o^yr(xi,x) : u : u/^+^(u, x, y; {u^-, x^, y^}; i) 

i 

-i^ Y^Vur j duj dxj dyL(xi, x)u/^+^(u, x, y; {u^, x^-, y^}; t) 

i 

As in the case of the one-point pdf no closed equation for the N-point pdf is obtained. The 
evolution equation for the N-point pdf contains the N-|-l-point pdf leading to a hierarchy of 
equations. 

In the dissipation term we can perform a partial integration which eliminates the formally 
defined operator L(xj,x). Additionally, one can recast the pressure and dissipation terms 
in a way which evidences Galilean invariance: 

^r({u,-,x,-,y,};t) +5:u, • VxJ^({u,-,x,-,y,};t) (28) 



= Z] • c?u dx <iyr(xi, x) : (u - Uj) : (u - Uj) 



/^+^(u,x,y;{uj,Xj-,yj};t) 

-'^J2'^ur J dujdxj rfy5(xi - x)Ax[u - Ui]/^+^(u,x,y; {uj-,Xj,yj};t) 

i 

Let us now introduce the notation 

({u,, X,, y,}; t) + Y: u. • V.J^({u,, X,, y,}; t) (29) 
= - Y^Vur J dn J dx J dyyl(xi-x,Ui-u)/^+^(u,x,y;{Uj,Xj,yj};t) 

i 

Here, a denotes the triple u, x, y. A is an operator which is related to the acceleration and 
is defined according to: 

A(xi - X, Uj - u) = -r(xi,x) : (u - Ui) : (u - Ui) + z/(5(xi - x)Ax[u - Uj] . (30) 

One may also add a random force field, which is Gaussian as well as 5-correlated in time. 
Then a diffusion term of the form 

^ E E Vu,Q(x„ X,) Vu,r ({a,}; t) (31) 

^ i j 

has to be included. 

It is important to formulate the following invariance properties of the hierarchy. Due 
to Galilean invariance /^({uj + c,Xj — ct,yj};t) solves the hierarchy provided that 
/^({uj, Xj, yj}; t) is a solution of the hierarchy. Neglecting the viscous terms, it can be shown 
that if /^({uj, Xj, yj}; t) is a solution, then also A^'^'^'^'^^/^dA'^u,;, A'^Xj, A'^y^}; At) solves the 
hierachy for each value of A, provided that 

5-7 = 1 . (32) 

The scale symmetry for an arbitrary value of S is obviously broken by the viscous term. 

As we have indicated above, the Eulerian N-point pdf can be calculated by the cor- 
responding Lagrangian pdf by integration with respect to the initial locations y^ of the 
particles. Integrating each equation of the hierarchy (^) we obtain a corresponding one for 
the Eulerian probability distribution function. This hierarchy has already been presented 
by Lundgren as well as Ulinich and Lyubimov 



Let us briefly comment on the question why a whole hierarchy of evolution equation arises. 
The mathematical treatment of fluid flow leads to a fleld theory. Considering only a flnite 
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number of fluid particles, therefore, yields a description with restricted information. This 
reduction of information shows up in the existence of an unclosed hierarchy of evolution 
equations for the joint velocity-position pdfs. A closed evolution equation can only be 
expected to arise when one approaches the continuum description, i.e. when one considers 
the full probability density functional defined in eq. (|TB|). 



IV. EVOLUTION EQUATION FOR THE PROBABILITY DENSITY FUNC- 
TIONAL 

This section is devoted to the derivation of a closed evolution equation for the probability 
density functional P[x(y), u(y), t]. Time differentiation yields: 

|p[x(y),u(y),t] = 

- I dy{< X(t,y) ■ ^-i-^i)[x(y) - X(t, y)]Z)[u(y) - U(t,y)] > 

+ < V{t, y) • -A-Z}[x(y) - X(t, y)]Z)[u(y) - U(t, y)] >} . (33) 
^u(y) 

Using the Lagrangian formulation of the Navier-Stokes equation ([T2|) we end up with the 
following relation: 

d ^ /• , _ 5 

iy) 



|p[x(y),u(y),t] + I rfyu(y) ■ ^P[x(y), u(y), t] 
= JdyJ dy'{^ ■ r(x(y),x(y')) : u(y') : u(y') 



(y) 

■^(x(y),x(y'))u(y')}P[x(y),u(y),t] . (34) 

We have arrived at a closed equation determining the evolution of the probability functional 
P[x(y), u(y), t]. The N-point probability functions /^({uj, Xj, yj}; t) can be obtained from 
P[x(y), u(y), t] by functional integration 

/^({u„x„y,};t) = J Du(y)Dx(y)5[ui - u(yi)]5[xi - x(yi)] 

(J[ujv - u(y7v)]5[xjv - x(yjv)]P[x(y),u(y),t] . (35) 



The hierarchy of evolution equations ( poD for the N-point pdfs is a projection of the functional 
equation ( pi]) onto the pdf of N fluid particles according to (|35D . We mention that, in 



principle, a projector formalism should be used to pass from the evolution equation for 
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the probability functional (|3^) to a closed equation for the projected N-point probability 
distribution. 



V. GENERALIZED FOKKER-PLANCK EQUATIONS 



If we investigate the problem of turbulence by a hierarchy of evolution equations such 



as (|29D we need to formulate suitable closure schemes. In the following we shall present 
a formulation of ( ^9]) which seems to be more suitable for that purpose. It is influenced 



by the so-called Lagrangian pdf method [0, [T^ , [jl5| , whose basic idea is to model the 
acceleration term by a stochastic force. Assuming Markovian properties the model pdf obeys 
a Fokker-Planck equation of the form 

d 



dt 



r({u,,x,,y,};f) + ^u, ■ VxJ^({u,,x„y,};t) (36) 

i 

= -EVu.-I^'({u,,x,,y,}) 

i 

+ Vu,I^'({u,-, X,, y,}) Vu,r ({u„ X,, y,}; t) . 



The question arises, what functional form of the drift term D^{{uj, Xj, yj}) and the diffusion 
matrix D^({uj, Xj, yj}) has to be chosen in order to obtain an accurate model of turbulent 



flows. Although the case of a single fluid particle seems to be well- investigated |jT3[, the case 



of several fluid particles has to be studied in more detail |T^, [|T8|, |T^. Interesting models 



for several fluid particles have been devised by Pumir et al. |T^, pO| 



Origionally, the random force method dates back to Oboukhov ||T6[, who suggested to 
use a Fokker-Planck equation without drift term and a constant diffusion term for the single 
particle case. 

In the following we shall derive a generalization of the Fokker-Planck equation directly 
from the hierarchy (pp|). To this end we consider equation to be a linear inhomogeneous 
equation, which can be solved in a straightforward manner (we consider the case of the 
(N+l)-point pdf): 



/^+^(u,x,y; {ui,Xi,y,};t) 



(37) 



- rt/t'e-(*-*')t"-^-+S»"-^''>l / da' 

Jtn J 
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A(x, - x', U,; - u') ■ Vu, + A(x - x', u - u') ■ Vu] 

i 

/^+^(u,x,y;u ,x',y'; {u^-, Xj, yj}; t') 



The first term stems from tlie initial condition. 

We remind tlie reader tliat we have considered the force free case. Nonrandom forces 
F(x, t) can be taken into account by a different evolution operator 



e 



-(t-t')u-Vx ^ yg-(i-i')u Vx-/"drF(x,r) Vu ^gg^ 



(Here, T denotes Dysons time ordering operator). Whereas the first evolution operator 
makes the replacement 

X ^ X - u(t - f) (39) 
the second evolution operator replaces 

X ^ X(x, t-f) , (40) 

where X(x, t — to) is the solution of the set of differential equations 

|^X(x, t-t') = U(x, t - t') 

|;U(x,t-t') = F(X(x,t-t'),0 (41) 

with the conditions 

X(x,0) = X 

U(x,0) = u . (42) 

Let us make some remarks on the external force. In three dimensional turbulence the 
force varies on the so-called integral scale, which is larger than the scales belonging to the 
inertial scale. That implies that during inertial time scales the relative motion of Lagrangian 
particles located within the inertial range is not influenced by the external force, i.e. by the 
mechanism how the turbulence is generated. Therefore, on time scales belonging to the 
inertial time scale the approximation ( PU] ) is sufficiently good for our purposes. 

Now we insert the obtained expression into the acceleration term of eq. (^). As a result 
we arrive at the hierarchy 

^r({u,,x„y,};t) + 5:u,-Vx,r(K,x„y,};f) (43) 
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= -Zl"^"- ■ / '^^'[D^(xi|{u,,Xi,yi};t,t')/^({ui,Xi,y,;};t')]x,=x,(x,,t-t') 

rt 

+ dt'[D^(xi,Xj|{ui,Xi,yi};t,t') ■ Vu,/^({ui, Xj, y^}; t')]x,=x,(x„i-t') 

- XI ■ / '^C^A(Xi - X, Uj - u) 

i 

/^+^(u,X(x,t - to),y; {ui,Xi(xi,t - to),yi};t = to) 

The last term stems from the initial condition. We want to point out that, in some sense, 
the present procedure is analogous to a projection operator formalism ||26|. It is well-known 
that the stochastic equations contain initial terms after projection. By some more or less 
sophisticated arguments, these initial terms are dropped. 

Our result ( ^3D takes the form of a generalized Fokker-Planck equation, where the gen- 
eralized drift term is given by 

D^{xi\{ni,Xi,yi};t,t')) = 

J da J dcr'A(x,; — X, Uj — u) 

M^j - x', - u') ■ Vu^ + A(i - x', u - u') ■ Vu] 

j 

p^+2(u',x',y';u,x,y|{uj,Xj-,yj};t')}x=x(x,t-t') • (44) 
The diffusion term takes the form 

D^{xi,±j\{ui,±i,yi};t,t') = 

J da J dcr'A(xj — x, Uj — u) A(xj — x', Uj — W) 
xp^+^(u', x', y'; u, x, y|{uj, Xj, y^}; t')li=x(x,t-t') 

A formally closed equation has been obtained by the introduction of the conditional proba- 
bility distribution 

r+2(a>;{a,};t) =p^+2(a';a|K};t)r({a,};t) , (45) 

{a denotes the triple u, x, y.) A successful description of Lagrangian turbulence statistics 
can be achieved if this conditional probability distribution can either be approximated or 
modeled in a suitable way. Thereby, the fundamental symmetries, i.e. Galilean and scale 
invariance (for u = 0), have to be retained. Furthermore, incompressibility of the fluid 
motion should be conserved. This requirement seems to be the major difficulty, since any 
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approximation has a consequence for the pressure term. However, only a correct treatment 
of the pressure term guarantees incompressibihty of the fluid. 



VI. SUMMARY 



We have formulated a hierachy of evolution equations for the Lagrangian N-point pdf 's 
in close analogy to the one for the Eulerian pdf's presented by Lundgren p3|] and Ulinich 



and Lyubimov [^. Due to the pressure and dissipative terms the N-point probability 
distributions couple to (N+l)-point distributions. The existence of a whole hierarchy of 
evolution equations is due to the fact that a field theory is described by a finite number of 
points. A closed statistical equation arises when one defines a probability functional. We 
have formulated an evolution equation for this functional which is the Lagrangian analog of 
Hopf's functional equation. 



Furthermore, we have tried to derive the Lagrangian pdf models ||T2[ , which are successful 
in modeling various aspects of turbulent flows by diffusion processes. Starting from the 
hierarchy of evolution equations for N-point pdfs we arrived at a generalized Fokker-Planck 
equation, i.e. a diffusion equation containing memory terms as well as a term steming from 
the initial condition. The generalized drift and diffusion coefficients are formally expressed 
by conditional probability distributions, so that the problem is not closed. In a following 



paper we shall adress the problem of formulating suitable closure approximations |25 



Acknowledgement: I gratefully acknowledge financial support from the Deutsche Forschungs- 
gemeinschaft within the project {Interdisziplindre Turbulenzinitiative) . Furthermore, I thank 
Joachim Peinke (Oldenburg) and Rainer Friedrich (Miinchen) for interesting discussions. 



[1] GI Taylor, Diffusion by continuous movement, Proc. London Math. Soc. Ser. 2, 20, 196 (1921) 
[2] L. F. Richardson, Atmospheric diffusion shown on a distance-neighbour graph, Proc. R. Soc. 

Lond. A 110, 709 (1926) 
[3] A. S. Monin, A. M. Yaglom, Statistical Fluid Mechanics: Mechanics of Turbulence, Volume 1 

(Cambridge, MA: MIT Press (1971)) 



14 



[4] A. S. Monin, A. M. Yaglom, Statistical Fluid Mechanics: Mechanics of Turbulence, Volume 2 
(Cambridge, MA: MIT Press (1975)) 



[5] 
[6] 

[7 
[9 

[lo: 
[11 

[12 

[13 
[14 

[15 
[16 

[17 

[18 

[19 

[2o; 



Prisch, Turbulence. The Legacy of Kolmogorov 

G. Falkovich, K. Gawedzki, M. Vergassola, Particles and Fields in Fluid Turbulence, Rev. 
Mod. Phys. 73, 797 (2001) 

A. La Porta, G. Voth, A. M. Crawford, J. Alexander, E. Bodenschatz, Fluid particle acceler- 
ation in fully developped turbulence. Nature 409, 1017 (2001) 

A. La Porta, G. Voth, F. Moisy, E. Bodenschatz, Using cavitation to measure statistics of 
low-pressure events in large Reynolds-number turbulence Physics of Fluids 12, 1485 (2000) 
Greg A. Voth, A. La Porta, A. M. Crawford, J. Alexander, E. Bodenschatz, Measurement of 
particle accelerations in fully developed turbulence, arXiv:physics/01 10027 vl 



N. Mordant, P. Metz, O. Michel, J.-F. Pinton, Measurement of Lagrangian velocity in fully 
developed turbulence, Phys. Rev. Lett. 87, 214501 (2001) 

S. Ott, J. Mann, An experimental investigation of the relative diffusion of particle pairs in 
three-dimensional turbulent flow, J. Fluid Mech. 402, 207 (2000) 

S.B. Pope, Lagrangian PDF methods for turbulent flows, Annual Review of Fluid Mechanics 
26, 23 (1994) 

S. B. Pope, Turbulent flows, Cambridge University Press (2000) 

5. B. Pope, Pdf methods for turbulent reactive flows. Prog. Energy Combust. Sci 11, 119 
(1985) 

B. Sawford, Turbulent relative dispersion, Annu. Rev. Fluid Mech. 33, 289 (2001) 

A. M. Obukhov, Description of turbulence in terms of Lagrangian variables. Adv. Geophys. 

6, 113 (1959) 

E. A. Novikov Two-particle description of turbulence, Markov property, and intermittency. 
Physics of Fluids A 1, 326 (1996) 

G. Pedrizetti, Quadratic Markov modeling for intermittent turbulence, Phys. of Fluids 11, 
1694 (1999) 

A. Pumir, B.I. Shraiman, and M. Chertkov, Geometry of Lagrangian Dispersion in Turbulence, 
Phys. Rev. Lett. 85, 5324 (2000) 

P. Castiglione, A. Pumir, Evolution of triangles in a two-dimensional turbulent flow, Phys. 
Rev. E 64, 056303 (2001) 

15 



[21] B.M.O. Heppe, Generalized Langevin equation for relative turbulent dispersion, J. Fluid Mech. 
357, 167 (1998) 

[22] E. Hopf, Statistical hydromechanics and functional calculus, J. Rat. Mech. Anal. 1 87 (1957) 
[23] T. S. Lundgren, Distribution functions in the statistical theory of turbulence, Phys. Fluids 
10, 969 (1969) 

[24] F. R. Ulinich, B. Ya. Lyubimov, Statistical theory of turbulence of an incompressible fluid at 

large Reynolds numbers, Zh. Exper. Teor. Fiz. 55, 951 (1969) 
[25] R. Friedrich, Conditional probability distributions for turbulence, preprint (2001) 
[26] R. Kubo, M. Toda, N. Hashitsume, Statistical Physics II (Springer- Verlag Berlin, Heidelberg, 

New York, Tokio, 1978) 



16 



